جغرافیا و مخاطرات محیطی. شماره سی و یکم. پاییز ۱۳۹۸ 


صص ۵۹۷ 
68 ۱ ۵ :101 


بر آورد میزان رواناب حوضه آبریز لنبران چای در استان آذربایجان شرقی: کاربرد مقایسه‌ای روش‌های 


محمدحسین رضائی مقدم" - استاد ژئومورفولوژی, دانشگاه تبریز؛ تبریز, ایران 
میراسدالله حجازی - دانشیار ژئومورفولوژی. دانشگاه تبریزه تبریز ایران 


عبدالله بهبودی - دانشجوی ژومورفولوژی. دانشگاه تبریز تبریز. ایران 


تاریخ دریافت: ۱۳۹۸/۶/۲۵ تاریخ تصویب: ۱۳۹۸/۹/۱۷ 


حکده 


شبیه‌سازی جریان رودخانه‌ها؛ پیش‌بینی رفتار هیدرولوژیکی حوضه‌های آبریز و درک درست از ملفه‌های 
مختلف چرخه هیدرولوژیکی برای برنامه‌ریزی و حفاظت از منابع آبی ضرورت دارد. از طرف دیگر نبود 
برآورد مناسب رواناب حاصل از بارش در حوضه‌های آبریز مدیریت بهینه منابع آب به‌ویژه مدیریت 
بهره‌برداری از سدها و شبکه آب‌رسانی را دچار مشکل می‌نماید. در اين پژوهش, بیلان آبی حوضه آبریز 
لنبران با استفاده از مدل ٩۷۷۸‏ شبیه‌سازی شد. برای واسنجی و اعتبارسنجی مدل ٩۷۷۸‏ از دو روش 
طلآنآت و 5/۳12 استفاده شد. نتایج حاصل با استفاده از چهار معیار ارزیابی مدل‌ها شامل معیار نش 
ساتکلیف» ضریب تبیین. نسبت باقیمانده میانگین مربعات خطا به انحراف از معیار داده‌های مشاهداتی و 
دو فاکتور 0-120100۲ و 1-10010۲ مقایسه شدند. نتایج حاصل از کاربرد این دو روش نشان داد که روش 
]907۳ دارای الگوریتم موترتری برای واسنجی و تعیین عدم قطعیت مدل در این حوضه است. بر این 
اساس با اطمینان بیشتری می‌توان از مدل ۹۷۷۸ واسنجی شده با الگوریتم 5]1۳1-2 در مدیریت منابع 
آب. کمی سازی سناریوهای تغییر اقلیم و تغییر کاربری در داخل حوضه استفاده کرد. 


۱ نویسنده مسئول: ۰۹۱۶۶۱۶۲۱۹۰ ۲ 011200ه) ۵ داعم 2ع۲ :اتمصرظ 


۰ جغرافیا و مخاطرات محبطی شماره سی و یکم 
۱- مقدمه 


حفاظت از منابع آب‌وخحاک ازجمله راهکارهای مهم مدیریت اصولی و صحیح حوضه‌های آبریز است؛ بنابراین 
فرسایش و نگرانی‌های مربوط به تغییر کاربری» اهمیت پیش‌بینی اثرات این اقدامات را بر منابع آب خاک بیشتر نمایان 
می‌کند. ازاین جهت ضرورت مدل‌سازی و ارزیابی نتایج آن‌ها ضروری می‌باشد. اغلب سازندهای تشکیل دهنده حوضه 
آبریز لنبران چای به علت ویژگی‌های خاص لیتولوژی دارای حساسیت بالایی نسبت به فرسایش بوده و عمدتاً از 
نفوذپذیری کمتری برخوردارند و باعث ایجاد رواناب سطحی می‌گردد. همچنین اختلاف ارتفاع زیاد حوضه. کمبود یا 
فقدان پوشش گیاهی. غیر اصولی بودن عملیات کشت تغییر مداوم کاربری اراضی؛ آتش زدن پوشش گیاهی. چرای 
بیش ازحد و لگدکوب شدن خاک توسط احشام اين امر را تشدید می‌کند؛ بنابراین برای مدیریت بهتر اطلاع از بیلان آبی 
حوضه‌های آبریز ضروری می‌باشد و به دلیل وقت‌گیر و پرهزینه بودن اندازه‌گیری اجزای بیلان آب» مدل‌سازی 
پدیده‌های هیدرولوژی در حوضه آبریز می‌تواند راه‌حل بهینه‌ای برای آن‌ها باشد (زارع کاریزی و طالبی» ۱۳۹۵)؛ بنابراین 
آگاهی دقیق از رفتار هیدرولوژیک حوضه‌های آبریز می‌تواند به شبیه‌سازی بهتر این محبط به‌منظور کنترل مولفه‌های 
اصلی بیلان آب کمک شایانی نماید (موذن‌زاده و همکاران ۱۳۹۵). از طرفی شبیه‌سازی فرآیندهای آب‌شناختی به دلیل 
هزینه‌های زیاد برای اندازه‌گیری و محدود بودن داده‌های در دسترس از اهمیت خحاص و ویژه‌ای برخوردار است ( 511۷۵۸ 
 21., 0‏ ونامحصتطگ ,2008 له اه ۲009002 ,2007 ,.له 61 وکاویان و همکاران, ۱۳۹۲)؛ اما کمبود آمان اطلاعات 
و منابع از موانع اصلی مسیر هر پژوهش می‌باشد بنابراین به روشی احتیاج است تا به‌وسیله آن بتوان آمار موخود را 
برای حوضه‌های آبریز بدون آمار يا مکان‌های تعمیم داد که اندازه‌گیری در آن‌ها امکان‌پذیر نیست ( .1 64 5هتقطوز 
0 مدل‌های زیادی برای شرح و پیش‌بینی هیدرولوژی آبریز پيشنهاد شده است که از نظر اهداف و مقیاس زمانی 
و مکانی بسیار متفاوت هستند (2010 ,81 61 56080). مدل‌های نیمه توزیعی هیدرولوژیکی مانند مدل ۹۱۷۸ اخیرا 
به‌طور گسترده‌ای برای مدیریت حوضه آبریز به کار گرفته می‌شوند (عینی و همکاران» ۱۳۹۷). توانایی مدل 5۷7۸۲ 
در شبیه‌سازی فرآیندهای هیدرولوژیکی پیچیده حوضه‌های آبریز در محیط 018 این مدل را نسبت به مدل‌های یکپارچه 
که در آن‌ها واحدهای کاربری بزرگ‌تر مبنای عمل هستند. متمایز ساخته است (اکبری؛ ۱۳۹۲). از میان داده‌های ورودی 
مدل 57۷۵۲؛ نقشه‌های توزیعی خاک و کاربری اراضی به دلیل در دسترس بودن. حساسیت کمتر مدل به تغییرات 
این نقشه‌ها و همچنین تغییرات نسبتاً کم آن در طول زمان و اساسا به دلیل پیشرفت‌های صورت گرفته در سنجش‌ازدور 
با اطمینان بالایی در کارهای مدل‌سازی به کار گرفته می‌شوند (2017 ,.21 61 01015500)). مدل ٩۲۷۸۲‏ یکی از مدل‌های 
است که تحقیقات گسترده‌ای در سطح جهان و ایران در زمینه‌های مختلف هیدرولوژی و مدیریت حوضه‌های آبریز 


صورت گرفته است که به برخی از کارهای انجام شده در ايران و جهان اشاره می‌شود. 
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(2015 ,۹600 تک مدا9) در پژوهشی با استفاده از مدل 81۸۲ نشان دادند که واسنجی یک مدل هیدرولوژیکی 
با استفاده از الگوریتم‌های متداول واسنجی خودکار لزوماً دقت مناسبی در پیش‌بینی متغیرهای هیدرولوژیکی در دوره 
اعتبارسنجی به وجود نخواهد آورد. (فرامرزی و همکاران ۲۰۰۹) با استفاده از مدل ۹۷۷۸۲ اجزای بیلان آب را در 
کل ايران شسبیه‌سازی کردند. سپس ارتباط آن‌ها را با تولید گندم با در نظر گرفتن عملکرد سدها و اقدامات آبیاری 
اراضی کشاورزی بررسی و مطالعه نمودند. (عباسپور و همکاران, ۲۰۰۷) برای شبیه‌سازی فرآیندهای موثر بر بیلان 
آب» رسوب و عناصر غذایی در حوضه آبریز رودخانه تور کشور سوئیس از مدل 7 استفاده کردند. نتایج نشان 
داد که شبیه‌سازی رواناب و نیترات بسیار حوب و شبیه‌سازی رسوب و فسفر نسبتاً وب صورت گرفته است. ( 1 
7 ,2 6۷) در تحقیق خود در حوضهای در غرب آفریقا با استفاده از مدل ٩۱۷۸۵‏ نشان دادند که تغییر کاربری 
حوضه‌ها از جنگل. مرتع و بوته زار به اراضی کشاورزی و يا مناطق شهری سبب تغییر واکنش هیدرولوژی حوضه 
می‌گردد. بیات و همکاران (۱۳۹۷) از مدل ٩۱۷۸‏ برای جذب داده - پیش‌بینی استفاده نمودند نتایج نشان می‌دهد 
که روند همزمان واسنجی - پیش‌بینی با اسستفاده از الگوریتم > برای این مدل باعث دقت پیش‌بینی می‌شود. 
گودرزی و فاتحی فر (۱۳۹۸) با استفاده از مدل ٩۷۷۸۲‏ در اثر تغییرات اقلیمی تحت سناریو 8.5 16۳ به پهنه‌بندی 
سیلاب در حوضه آذرشهرچای پرداختند. شفیعی مطلق و همکاران (۱۳۹۷) تأثیر کاربری اراضی را بر رواناب رودخحانه 
مارون با استفاده از مدل ٩۷۷۸۲‏ بررسی نمودند نتایج نشان داد که در منطقه سناریوی بدبینانه اتفاق است و با ادامه 
روند تخریبی مقادیر مختلف رواناب افزایش و نفوذپذیری و آب گذری به آبخوان‌های سسطحی و عمیق کاهش 
می‌یابد. بصیری و همکاران (۱۳۹۳) پژوهشی با عنوان تأثیر مقیاس نقشه خاک بر مدل ٩۱۷۸۲‏ در شبیه‌سازی روند 
جریان رواناب ماهانه در حوضه آبریز تویسرکان انجام دادند. نتایج نشان داد که تنوع خاک‌های منطقه مطالعاتی نقش 
مهمی در نتایج مدل ایفا می‌کند. 

(عارفی اصل و همکاران. ۱۳۹۲؛ بسالت پور و همکاران. ۱۳۹۳؛ نظری پویا و همکاران۱۳۹۶؛ ناصرآبادی و 
همکاران, ۱۳۹۵؛ غلامی و همکاران: ۱۳۹۲: محمدی و همکاران, ۱۳۹۲؛ اعلمی و همکاران ۱۳۹۷؛ حبیبی و همکاران, 
۳۷ در حوضههای متفاوت در نقاط مختلف ايران با مدل ٩۱۷۸۲‏ به شبیه‌سازی رواناب پرداختند. به‌طورکلی 
هدف از انجام این پژوهش‌ها عمدتاً برای برآورد میزان رواناب و رسوب جهت بررسی تغییرات کاربری و مدیریت 


بهینه منابع آبی در حوضه آبریز بوده است. 
۲- مواد و روش 
۱1-۲- معرفی منطقه مورد مطالعه 


زیر حوضه لنبران در سمت ضلع شرقی رودخانه اهرچای واقع شده است. مساحت این زیر حوضه ۲۰۱۱۸ 


هکتار می‌باشد و سرشاخه اصلی رودخانه اهرجای محسوب می‌شود. محدوده مورد مطالعه در بین مختصات 


1 جغرافیا و مخاطرات محیطی ی 
جغرافیایی ۱۷ ۲۰2 47۳ تا 4۵ [۳۵ "۶7 طول‌های شرقی و 1٩‏ ۲۰2 ۳۸ ۲۳2 ۳۵ ۳۸ عرض‌های شمالی 


استان آفربایجان شرقی واقع شده است (شکل ۱). 


00000 30000 560000 


00000 
0000 


حو ضه آب ری ز ارس 


0000 
00 


حو ضه آبٍ ری ز اهر 


000 
00 


۱ 


00000 30000 560000 


دوتات 
42,0 0 21,00 


شکل شماره ۱- موقعیت منطقه مورد مطالعه. 

۲-۲-مدل 5۷۷۸۲ 

71 نمونه‌ای از مدل‌های فیزیک پایه است که با حل معادلات اساسی فیزیک به شبیه‌سازی فرآیندهای سیستم 
آبریز می‌پردازد. این مدل به لحاظ مقیاس فضایی نیمه توزیعی و از نظر زمانی مدلی پیوسته است. کوچک‌ترین واحد 
کاری در این مدل واحد واکنش هیدرولوزیک 11817 است که از ترکیب نقشه‌های طبقات شیب خاک و کاربری 
اراضی حاصل می‌شود. آب موجود در خاک رواناب سطحی, رسوب و عناصر شیمایی ابتدا برای هر 11181 و سپس 
هر زیر حوضه و در آحر کل حوضه آبریز محاسبه می‌شود. شبیه‌سازی سیستم آبریز در مدل ٩۷۷۸۲‏ را می‌توان به دو 
بخش کلی تقسیم نمود: فاز زمینی و فاز آبی. فاز زمینی مربوط به فرآیندهای سسطح زمین و ورود آب» رسوب و 
عناصر شیمیایی به آبراهه اصلی هر زیر حوضه است. فاز آبی فرآیندهای آبراهه‌ها و کانال‌های جریان شامل حرکت 
آب» رسوب و مواد شیمیایی را شبیه‌سازی می‌کند (2005 ,21 6 100). شبیه‌سازی بخش زمینی چرخه 
هیدرولوژیک در مدل 5۷۸۲ بر پایه رابطه بیلان آبی رابطه (۱) صورت می‌گیرد. 


تلا معمممومگ ۲۱۷۲۵۱0212 1 
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رابطه (۱) ار 9 - مووو۷ - ی - ار ی وک < 5۱۷ 
که در این معادله ٩۷۷۲‏ مقدار نهایی رطوبت خاک (میلی‌متر)؛ 9۷۷0 مقدار اولیه رطوبت خاک (میلی‌متر» ] زمان 
(روز» وکا مقدار بارش در روز 1 ام (میلی‌متر)؛ ۶«) مقدار رواناب در روز 1 ام (میلی‌متر» مت مقدار تبخیر و تعرق 
در روز 1 ام (میلی‌متر)؛ «معه ۷۷ مقدارنشت که در روز 1 ام از نیمرخ خاک وارد منطقه غیر اشباع می‌شود (میلی‌متر)؛ 
«ع) مقدار آب بازگشتی در روز 1 ام (میلی‌متر) می‌باشد. در مدل ٩۷۷۸‏ پروفیل خاک به چندین لایه تقسیم 
می‌شود و فرآیندهای حرکت آب در خاک شامل رواناب سطحیی نفو تبخین جذب توسط گیاهان. جریان جانبی و 
نفوذ عمقی به سفره‌های آبی کم عمق و عمیق می‌باشد. شبیه‌سازی رواناب سطحی با روش اصلاح شده 902 (در 
صورت وجود بارش روزانه) يا روش گرین - آمپت " (در صورت وجود داده‌های بارش ساعتی) انجام می‌شود. برای 
شبیه‌سازی حرکت آب در داخل لایه‌های خاک از روش روندیابی مخزن استفاده می‌شود. جریان اشباع ب‌صورت 
مستقیم و جریان غیر اشباع بین لایه‌های خاک به‌طور غيرمستقيم با محاسبه جذب گیاهی و تبخیر از مسطح خاک 
شبیه‌سازی می‌شود. جریان آب به سمت لایه‌های پایین زمانی اتفاق می‌افتد که میزان آب در یک لایه از خاک از 
ظرفیت مزرعه تجاوز کند و لایه زیرین نیز غیراشباع باشد. جریان جانبی با استفاده از روش روندیابی مخزن جنبشی 
(1983 .ة 64 50180). بر مبنای درجه شیب طول شیب و هدایت هیدرولیکی اشباع, مدل‌سازی می‌شود ما بقی آب 
نفوذ يافته, به سفره‌های کم عمق زیر زمینی وارد می‌شود. برای برآورد تبخیر و تعرق پتانسیل» روش‌های هارگریوز - 
سامانی » پریستلی - تیلور * و پنمن - مانتیث * موجود می‌باشد. مدل 5۷۷۸۲ تبخیر از سطح خاک و تبخیر از گیاه را 
ب‌صورت مجزا مطابق روش ربچی (1972 ,8696 محاسبه می‌کند. تبخیر از سطح خاک با یک تابع نمایی از عمق 
خاک و میزان آب خاک برآورد می‌شسود. تبخبر از گياهان نیز با یک تابع حطی از تبخبر و تعرق پتانسیل شساخص 
سطح برگ عمق ریشه و میزان آب خاک مدل‌سازی می‌شود. میزان رطوبت لایه‌های مختلف خاک در هر زمان و 
مکان با در نظر گرفتن مجموع عوامل فوق برآورد می‌شسود. پس از تعیین رواناب سسطحی با روش شماره منحنی یا 
گرین - آمپت. مقدار رواناب وارد شده به کانال آبراهه محاسبه می‌شود. معادله مانینگ" برای تعیین سرعت جریان و 
روش ذخیره متغیر یا ماسکینگام " برای روندیابی جریان آب در کانال استفاده می‌شود. در ضمن, تلفات انتقال آب در 
مسیر کانال و تأخیر زمانی رواناب سطحی برای حوضه‌های آبریز بزرگ که زمان تمرکز آن‌ها بیش از یک روز است 
نیز در نظر گرفته می‌شود. 
68 60856۳۷2110۴8 9011 1 
0۲6۵8۵۸ 2 
تصحصصصصح 65-5 ۳۱2۲۵۲6۵2۷ 3 
۲ -۳۲۱۵۹01۷ 4 
۳۵۴۸۳۱۵۸-1 5 


م10 عمتصصه]۱۷ 6 
دنام م۱۷۲۹ 7 
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۱-۲-۲- روش تعیین عدم قطعیت ,۶1,۳ 

این روش اولین بار توسط (.,1992 ,13210 تک 136۷67) به‌عنوان روشی جدید که معمولاً با مدل‌های شبیه‌سازی 
زیست‌محیطی به کار می‌رود. معرفی گردید. در این روش عدم قطعیت شامل سه نوع عدم قطعیت داده‌های ورودی» 
ساختار مدل و عوامل خواهد بود زیرا مقدار تابع درست نمایی با مجموعه عوامل در ارتباط بوده و به‌طور ضمنی 
منشأ همه حطاها و تأثیر کواریانس مقادیر عوامل بر عملکرد مدل را منعکس می‌کند. روش 7ارآ) شامل سه مرحله 

۱- بعد از تعریف تابع درست نمای تعمیم یافته (1)0 تعداد زیادی مجموعه عوامل به‌صورت تصادفی از توزیع 
احتمال پیشین انتخاب می‌شود و بر اساس مقدار تابع درست نمایی و معبار حدی تعریف شده هر مجموعه از عوامل 
به‌صورت هم رفتار با غیر هم رفتار طبقه‌بندی می‌شود. 

1,0 

لا 

(۵) یر 


رابطه (۲) 


داده می‌شود. 

۲-۲-۲- روش تعیین عدم قطعیت 5/۳1-2 

در روش 51/۳12 همانند 01/1710 سنجش میزان عدم قطعیت در مدل‌سازی شامل عدم قطعیت در ورودی‌هاء 
مدل مفهومی عوامل می‌باشد و توسط معیار 100107- ۲ ارزیایی می‌گردد که بیانگر درصد داده‌های اندازه‌گیری شده 
که درون باند عدم قطعیت ۹۵ درصد قرار گرفته‌اند می‌باشد. مراحل انجام کار به‌صورت زیر می‌باشد: 

۱- در مرحله اول تابع هدف (8)0 و بازه عوامل سم .من«حه0 ) تعریف می‌شوند. 

۲- از روش نمونه برداری لأتين مکعبی در بازه تعریف شده در مرحله قبل. جهت تولید مجموعه عوامل و مقدار 
تابع هدف به ازای هر مجموعه از عوامل استفاده می‌شود. مقادیر تابع هدف با استفاده از ماتریس ژاکوبین و ماتریس 
و کواربانس عوامل 6 مطابق رابطه زیر ارزیابی می‌گردد ("ی واریانس مقادیر تابع هدف از 130 بار اجرای مدل و 9 
تعداد عوامل). 


۸ 
رابطه () ی ی که زیربت کب .رز 


٩۵-۳‏ درصد از بازه پیش گوبی شده هر کدام از عوامل از رابطه زیر محاسبه مور کر قگ: 


سال هشتم برآورد میزان رواناب حوضه آبریز لنبران چای در استان ... 10 


عی< 
اطه (۶ س سس 
راب 2 77 0025 2 بر باب 7- ره 
7+ وم 
رل 0۵25 0۵ بو ۳2 07272607 زر 


که در آن 6 مقدار پارامتر (6 در تابع هدف بهینه می‌باشد ۷ درجه آزادی یا 11۳0 می‌باشد. 
۳- بازه 9300 محاسبه شده سپس دو شاخحص 0-120107 درصد داده‌های مشاهداتی که در بازه 950۷ قرار 
می‌گیرند و شاخص "100007 1 از رابطه زیر محاسبه می‌شوند: 
و و 7 1 
‌ 2 فك 946 0( 1 2« ۳۹ 


وورط) 
که در آن »دوز و ۳2596 معرف مرز بالا و پایین عدم قطعیت 95000 و «:0) انحراف از معیار داده‌های 


[اطه( 9 < 0610۳ - 7 


مشاهداتی می‌باشد. بهترین نتیجه زمانی خواهد بود که ۰ درصد داده‌های اندازه‌گیری با 9500۷ پوشانده شود با 
بعبارت دیگر میزان ۲ برابر ۱ شود و 1-10010 کمتر و به صفر نزدیک باشد. در صورتی که مقادیر دو شاعص 
مذکور رضایت بخش باشد. بازه تعریف شده عوامل به‌عنوان توزیع احتمال پسین خواهد بود در غیر این صورت بازه 


جدید عوامل بر اساس رابطه زیر محاسبه شده و تمام مراحل ۱ تا ۶ تکرار خواهد شد. 


زانط( 
را 9 ون 0 - 0 
0 ز حقحط. [ صتص . ز 7 #- 
۰ و 7 :5 وس ی 9 
ری 9 ین 09 - 
راتس لت بت ری پرو ور یش 6 ی بای ۳ 


در روش ٩1/۳712‏ فرض می‌شود که هر پارامتر ناشناخته به‌طور یکنواخت در یک دامنه با عدم قطعیت معین 
توزیع شده است لازم به ذکر است در هر دو روش. درجه عدم قطعیت از دو فاکتور 0-126107 و 106107 -1 محاسبه 
می‌گردد. 

۳-۲- معیارهای عملکرد 

جهت تعیین میزان کارایی دو روش از ۶ معیار عملکرد مدل‌ها استفاده شده است. مقدار 187 بیانگر ارتباط خحطی 
بین داده‌های شبیه‌سازی شده و واقعی بوده که مقدار آن بین صفر تا یک متغییر است. ضریب 8 بزرگی نسبی 
واریانس باقیمانده را در مقایسه با واریانس داده‌های مشاهداتی تعیین می‌کند. (1999 ,له 61 ۷1000۷110۷) در بررسی 
تعیین عملکرد مدل‌ها به این نتیجه رسیدند که اگر مقدار ضریب ٩‏ بیش از ۰/۷۵ باشد شبیه‌سازی خوب و مقادیر 
بین ۰/۳۹ و ۰/۷۵ شببه‌سازی قابل قبول و کمتر از ۰/۳۱ غبر قابل قبول خواهد بود. پارامتر 9۴ با استاندارد نمودن 
مقادیر باقیمانده میانگین مربعات خطا ۷15۳[ با استفاده از انحراف معیار مقادیر مشاهداتی به دست ی ی برای این 
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معیار مقادیر بالای ۰/۷ قابل قبول می‌باشد. معیار 918 که در آن ضصریب تعیین 12 در ضریب خط رگرسیون بین 
داده‌های مشاهداتی و شبیه‌سازی شده ضرب شده است. مقدار پهینه آن یک می‌باشد با افزایش گام زمانی معیارهای 
ارزیابی سخت‌تر خواهد شد و بالعکس (2007 ,2 )6 [ک1۷1013). بر این اساس با توجه به اینکه مدل در گام زمانی 
روزانه اجرا شده است معیارهای مذکور می‌تواند با سخت گیری کمتری نسبت به گام زمانی ماهانه بررسی شود. 
۲-- داده‌ها 
داده‌های مورد استفاده در این پژوهش شامل نقشه توپوگرافی ۱:۲۵۰۰۰ سازمان نقشه برداری» کاربری اراضی» 
خاک بارش, دمای حداقل و حداکثر سرعت باد. تشعشع خورشیدی و رطوبت نسبی روزانه ایستگاه‌های سینویتیک 


اه تبخیر سنجی ورزقان و باران سنجی اویلق و دبی ایستگاه کاسین می‌باشد که (جدول ۱). 


حدول ۱- داده‌های مورد استفاده در مدل ٩۲۷۸۲‏ و واسنحی مدل. 


ردیف داده منابع بازه اطلاعاتی توضیحات کاربرد 
۱ مدل رقومی ارتفاعی نقشه برداری وود ۱9۰ ورودی 5۷۷۸ 
۲ نقشه کاربری اراضی منابع طبیعی ود تصاویر لندست ورودی 5۷۷۸۲ 
۳ نقشه خاک منابع طبیعی ود ۱ ورودی 5۷۷۸۲ 
‌ داده‌های هواشناسی هواشناسی ۷۰۹ ایستگاه اهر و ورزقان | ورودی 5۲7۸۲ 
۵ داده‌های بارندگی سازمان آب و هواشناسی ۹-- ۲۰۰۰ اویلق ورودی ۷۷۸۲ 
1 دبی ماهانه سازمان آب ۲۰۰۰۹ ایستگاه کاسین واسنجی ]5۷۷۸ 


در حوضه آبریز لنبران چای از نقشه خاک و نقشه کاربری اراضی تفضیلی تهیه شده توسط منابع طبیعی استان 
آذربایجان شرقی استفاده شده است براساس این نقشه‌ها در کل حوضه هفت کلاس بافت خاک قابل تفکیک می‌باشد. 
۷ زیر حوضه و 1٩‏ واحد هیدرولوژیکی تقسیم گردید. از سال آماری ۰ تا ۲۰۰۳ برای گرم کردن مدل از ۷۰ 


الم برای واسنجی و ۶ سال برای اعتبار سنجی استفاده شده است. 


۵- نتایج و بحث 
مقر ازیو هه شیم تاش یک وه اک هه ۳ تا 
اساس گام زمانی ماهانه صورت گرفت. به‌منظور تعیین درجه حساسیت پارامترهای جریان در مدل 5۷۷۸۲ با استفاده 


از دو روش 2 و تاأنابآت) حساسیت سنجی برای ۳۲ پارامتر انتخابی» انجام گردید. با استفاده از نتایج صحت 
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سنجی در مورد حذف پارامترهایی که دارای درجه حساسیت کمتر می‌باشند از فرآیند واسنجی تصمیم گیری می‌شود 
که درنهایت پارامترهای با درجه حساسیت کمتر حذف گردید و ۱۵ پارامتر که حساسیت بیشتری داشتند انتخاب 
گردید. بین ۱۵ پارامتر حساس مدل ۹۷۷۸۲ به عوامل چکالی توده حاک. هدایت هیدرولیکی خاک در حالت اشباع 
ضریب تبخیر آب زیر زمینی» ثابت تخلیه جریان از کانال؛ درجه حرارت ریزش برف درجه حرارت ذوب برف 
حساسیت بیشتری دارد. مهم‌ترین عواملی که برای واسنجی و تعیین عدم قطعیت مدل استفاده شده است به همرا 
مقادیر اولیه و مقدار واسنجی شله در (جدول ۲) آورده شده است. نتایج نهایی معیارهای ارزیابی مدل در دوره 
واسنجی و اعتبار سنجی رواناب در (جدول ۲ ارائه شده است. طبق نظر (2007 ,له اه لعقت:0۱0 ( ,.ل2 اه ۸۵۱616 
27 ) کرمی و خطیبی (۱۳۹۸) اگر ضریب کارایی مدل نش - ساتکلیف در ارزیابی مدل بین ۰/۷۵- ۰/2۵ باشد 
نتیجه شبیه‌سازی خوب است (جدول ۶). همچنین با مقایسه ضریب تبیین بین رواناب واقعی و شبیه‌سازی شده در 
دوره واسنجی و اعتبار سنجی می‌توان گفت رابطه معنی دار و خوبی بین آن‌ها وجود دارد که این موضوع با نتایج 
رضا زاده و همکاران (۱۳۹۷» عینی و همکاران (۱۳۹۷) و محمدی و همکاران (۱۳۹7) مطابقت دارد. هر چند کاربرد 
روش ٩11۳12‏ آسان می‌باشد ولی کاملاً خود کار نیست در استفاده از این روش لازم است توزیع پسین عوامل در هر 
تکرار جهت اصلاح بازه پیشنهادی جدید (توزیع پیشین) تکرار بعهدی بررسی گردد. واسنجی مدل برای رواناب در 
یک مرحله با تعداد ٩۰۰‏ شبیه‌سازی در سه تکرار انجام گردید فرآیند واسنجی زمانی به پایان می‌رسد که براساس تابع 
هدف ضرایب مورد نیاز جهت ارزیابی در حد قابل قبول باشند (جدول ۳. اشکال ۲ و ۳ دبی شبیه‌سازی شده و 
مشاهده‌ای را به ترتیب در واسنجی و اعتبارسنجی نشان می‌دهد. تمام معیارهای ارزیابی مدل در شبیه‌سازی رواناب 


در حد مجاز می‌باشد. 


حدول ۲- پارامترهای حساس به همراه توزیع پیشین و پسین و مقدار بهینه آن‌ها جهت استفاده در واسنحی 


رواناب در دو روش 8۳1-2 و 61,71 در حوضه آبریز لنبران. 


5912 طنااه 
پارامتر علامت اختصاری 
حداقل | حداکثر | بهینه حداقل | حداکثر | بهینه 
شماره منحنی 50 ۳002 ۰/۲ ۲ | ۱/۱۸ ۵ ۲۵ ۹- 
عمق آستانه شروع جریان زیر زمینی 000۲ ۷ ۲ ۰ ۱/۶ 15۰ ۰۰ ] ۱۰۸۰ 
ضریب انتقال از سفره کم عمق به عمیق 60 ۱۷ ۱ ۰ ۲ ۰/۲ ۱ ۹ 
ثابت تخلیه جریان از کانال هر( ]۷ ۱ ۰ ۰/۵ ۰ ۰/۵ ۰/۱ 
فاکتور جبران تبخیر از خاک ۳۹00 ۷ ۱ ۰/۱ ۰/۸۲ ۰/۲ ۱ ۳۳ 
هدایت هیدرولیکی خاک در حالت اشباع 1707/1 (101 ,50 ۲ ۰/۸ ۰/۸ ۵- ۵ ۲۵ ۳ 
متوسط شیب واحدهای هیدرولوژیک نصا ۲ ۰/۲ ۱ ۰/۱ 0 ۳۹ ۹/۲ 


پارامتر 


ظرفیت آب قابل دسترس در لایه خاک/ 11110 
چگالی توده اک 
درجه حرارت ریزش برف/ سانتی گراد 
ضریب ذوب برف /سانتی گراد 
فاکتور جریان برداشت آب گیاه 
ضریب تبخیر آب زیر زمینی 
ضریب زبری دأمنه 


درجه حرارت ذوب برف/سانتی گراد 


ادامه حدول ۳۲ 


علامت اختصاری 


۵ ,0 ۲ 
( ال ۲ 
۳ ۱۷ 
۲ ۲۷ 

۱۷ ۳۵0 
۷0۱۷۷ ۷ 
۷ 0۷ 
۲۷ ۳ 


بهینه حداقل 

-۲۵ [ ۳ 
0 ۰/۵ 
-ِ 1/۸ 
۰ ۲/۸ 
۰/۸ / 
۰ ۰/۸ 
۰/۲ ۰/۹۸ 
۰ 1/۱ 


ور 0 نی 

حداکثر | بهینه 
۳۵ ۰/۱ 
را ۳ 
۲ ۰/۹ 
۷۰ 5/۳ 
۱ ۰/۳۱ 
۱/۵ ۱/۷ 
۱ ۰/6۵ 
۳ ۳/۸ 


جدول ۳- نتایج نهایی معیارهای ارزیابی مدل در دوره واسنجی و اعتبار سنجی رواناب در حوضه آبریز لنبران. 


۳۲ 
۰/۸۱ 
۳۹ 
۰/۹۲ 


۱/۵ 


دوره آماری 
۰ سال 
۶ سال 
۰ سال 


۶ سال 


حدول ع- مقادیر و طبقات شاخص‌های 12 و د۷. 


دا 1 ۱ 
۳ ۰/۸۹ ۸۳ 
۰/۷۸ ۷ ۰/۹۱ 
۰/۷ ۳/۵ ۰/۸ 
۰/۲ ۰/2۹ ۰/۹۱ 

طبقات 
ناکار امد 
رضایتبخش 
خوب 
خیلی 
ناکار آمد 
قابل قبول 
رضایت 
وب 
خیلی 


۱/۵ 


۰/۵ -۵ 


۰/۱۵ -۵ 


۱/۸۷۵ -۱ 


منبع: موریاسی و همکاران (۲۰۰۷) آی آل (۲۰۱۷) و کرمی و بیاتی خطیبی (۱۳۹۸). 


الگوریتم 


۱0۳2 


نان 


شاخص 


۳ 


5 
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هر چند کاربرد روش 51777172 آسان می‌باشد ولی کاملاً حود کار نیست در استفاده از این روش لازم است توزیع 
پسین عوامل در هر تکرار جهت اصلاح بازه پیشنهادی جدید (توزیع پیشین) تکرار بعدی بررسی گردد. واسنجی مدل 
برای رواناب در یک مرحله با تعداد ٩۰۰‏ شبیه‌سازی در سه تکرار انجام گردید فرآیند واسنجی زمانی به پایان می‌رسد 
که براساس تابع هدف ضرایب مورد نیاز جهت ارزیابی در حد قابل قبول باشند (جدول ۳). اشکال ۲ و ۳ دبی 
شبیه‌سازی شده و مشاهده‌ای را به ترتیب در واسنجی و اعتبارسنجی نشان می‌دهد. تمام معیارهای ارزیابی مدل در 


شبیه‌سازی رواناب در حد مجاز می‌باشد. 


16 
14 
٩2 -‏ 
1۷40 
3 
9 
4 
2 
_ ی 30 یس هت 1 مستت 3 ات۳ 3 زو ده پپبه-۵ 
6٩ 686 71 76 81 86 81 86 161 108 111 8‏ 56 51 48 41 36 31 26 21 16 1 6 1 
ماهها 
[اک هد متسه سوتزی مایم یی 


شکل ۳- نمودار دیی واقعی و شبیه‌سازی شده به روش 2 در مرحله واسنجی. 


دبی | متر مکعب 
5 ۵ ه + ه و ج ب نم د ه 


جط 


مس و 
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شکل ۶- نمودار دیی واقعی و شبیه‌سازی شده به روش 2 در مرحله اعتبار سنحی. 


۷۰ جغرافیا و مخاطرات محبطی شماره سی و یکم 


تابع هدف ولا برابر / انتخاب ان این روش با تعداد نمونه گیری 3 ۰ و ۱۹۷۰۰ انجام گردید نتایج 
نشان داد با تعداد نمونه گیری ۸۰۰۰ عملکرد بهتری نشان می‌دهد. اشکال ۶ و ۵ دبی شبیه‌سازی شده و مشاهده‌ای 
را به ترتیب در واسنجی و اعتبارسنجی نشان می‌دهد. تمام معیارهای ارزیابی مدل در شبیه‌سازی رواناب در حد مجاز 


می‌باشد ولی مقایسه دو روش نشان می‌دهد روش ٩10712‏ نسبت روش تالآ[ در این حوضه عملکرد بهتری دارد. 


آدبی شبیه ساز ی شدقات دبی ر العی‌ رس 


شکل ۵- نمودار دبی واقعی و شبیه‌سازی شده به روش ناه در مرحله واسنجی. 
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شکل 7- نمودار دبی واقعی و شبیه‌سازی شده به روش ارآ در مرحله اعتبار سنحی. 


مقایسه گرافیگی, مقادیر دبی واقعی و شبیه‌سازی شده در مقیاس ماهانه نان می‌دهد که مدل ۹۷۷۸۲ در 


شییه‌سازی وواناب» جربان‌های حلاگر را کر از سل واقعی برآوود می کند. اکتر جریان‌های که مدل قاذر به 


شبیه‌سازی آن‌ها نبوده در اواخر زمستان و فصل بهار بوده‌اند یکی از دلایل ضعف مدل در برآورد حداکثر جریان 
رونت ام دهاز ملالن سازمان حتاطی عفا ی ور ماس به روانات اس که این مدل در روازان‌ساصا اروت 
ترف زا شسبیاسازی ی کنط اقران وسکاران (۱۳۸) ی اتسار و هسکاوان:(۱۳۹۱) گرازش کرفاند که اک 
تخریان‌هایی کهرمللن قادر بهاشحیستازی آن‌ها نبزدهفر فصتا بهار اعاق:می‌افتز مت 5۳/۸۲ پارقن زا با انس‌فاده از 
متوسط دمای روزانه به‌صورت باران يا برف تقسیم بندی می‌کند. با توجه به ضعف مدل در شبیه‌سازی رواناب ماهانه 
اسفند فروردین و اردیبهشت و تقدم جریان شبیه‌سازی شده فروردین در بعضی از سال‌ها می‌توان نتیجه گرفت که 
مدل قادر به شبیه‌سازی کامل ذوب برف نبوده و مقادیر حداکثر رواناب را برای این حوضه کوهستانی در ماه‌های 
مذکور شبیه‌سازی نکرده است. در این مدل هنگامی که رطوبت موجود در هر لایه خاک از گنجایش زراعی بیشتر 
باشد مقداری از این آب اضافی باعث ایجاد جریان جانبی می‌شود و قسمتی هم به لایه بعدی نفوذ می‌کند؛ اما هنگامی 
که لایه‌ای از خاک یخ زده بافد مدل فرض می‌کند که هیچ جریان جانبی اتفاق نمی‌افتد و لایه حاک می‌تواند آب 
اض‌افه بر گنجایش زراعی را تا زمانی که اشسباع شود در خود نگه دارد. سپس فرض می‌کند که رطوبت مازاد بر 
گنجایش زراعی می‌تواند از لایه اسباع یخ زده به لایه بعدی نفوذ کند. (2004 .,500۵:001607 ۵20 0[902) بیان نمودند 
که مدل 5۱۷۸۲ برای شبیه‌سازی وقایع شدید طراحی نشده و معمولاً حداکثر جریان را کمتر از مقدار واقعی برآورد 
می‌کند (2004 ,ن0متهطمنط؟ ک تط)) در پژوهشی در مریلند نشان دادند که مدل 51۷۸ قادر نیست شرایط 
هیدرولوژیک را در دوره‌های تر سالی شبیه‌سازی کند. ازجمله علل دیگر ضعف مدل در نظر نگرفتن درصد تراکم 
پوشش گیاهی است که در اين مدل صرفاً به نوع پوشش گیاهی اشاره شده است. از دلایل دیگر ضعف مدل در 
شبیه‌سازی جریان در بعضی از ماهها ی سال می‌توان به تعداد داد‌های کم جریان, و فقدان داده‌های ورودی لازم برای 


شبیه‌سازی آب زیر زمینی اشاره کرد. 


1- جمع بندی 

با توجه به نتایج حاصل از حساسیت سنجیء چگالی توده حاک, هدایت هیدرولیکی خاک در حالت اشباع ضریب 
تبخیر آب زیر زمینی, ثابت تخلیه جریان از کانال, درجه حرارت ریزش برف» درجه حرارت ذوب برف حساسیت 
بیشتری دارد. هم چنین بررسی نتایج حاصل از کاربرد دو روش 5117712 و 01.17 در حوضه آبریز لنبران نشان داد 
که بر اساس معیارهای ارزیابی ضریب تبیین و نشساتکلیف. روش ٩10۳12‏ هم در مرحله واسنجی و هم در مرحله 
اعتبار سنجی نتایج خوبی نسبت به روش تاأنآنآ) در این حوضه دارد. هم چنین هر دو روش نمی‌توانند دبی‌های اوج 
را بخوبی پیش‌بینی کنند زیرا هر دو روش در دبی‌های اوج دبی. را نسبت به دبی واقعی کمتر پیش‌بینی کرده‌اند. از 


دلایل ضعف مدل در شبیه‌سازی جریان در بعضصی از ماه‌های سا می‌توان به شیه‌سازی ضعیف ذوب برف» ناسا زگاری 
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۷۲ جغرافیا و مخاطرات محبطی شماره سی و یکم 


فرضیات مدل در انتقال جریان در لایه‌های یخ زده و اشباع تعداد کم داده‌های جریان و فقدان داده‌های ورودی لازم 
برای شبیه‌سازی آب زیر زمینی اشاره نمود. هر چند کاربرد روش ]نآ آن) در مقایسه با روش 511۳12 آسان‌تر می‌باشد 
ولی در روش 01.1۳8 تعداد شبیه‌سازی جهت رسیدن به جواب بیشتر است؛ بنابراین با توجه به نتایج بدست آمده 
درنهایت می‌توان گفت که روش 9171772 با تعداد شبیه‌سازی کم در مقایسه با روش نانآ در این حوضه نتایج 
بهتری نشان می‌دهد. برای مشحضی شدن کارایی بهتر مدل. پيشنهاد می‌شود اولاً از آماز ایستگاه‌های دبی سنجی متعند 
و دبی روزانه به جای دبی ماهانه» در صورت وجود در یک حوضه استفاده شود تا نتایج آن‌ها با هم مقایسه شوند. ثنیً 
برای رسیدن به نتایج مطلوب این مدل بصورت مقایسه‌ای با سایر مدل‌های شبیه‌سازی در این حوضه و حوضه‌های 


همجوار استفاده گردد. 
کتابنامه 


اخوان. سمیرا» عابدی کوپایی. جهانگیر. موسوی, سید فرهاد. عباس پور کریم. افیونی. مجید و اسلامیان» سید سعید؛ 
۰ تخمین «آب آبی» و «آب سبز» با استفاده از مدل ٩۷۷۸‏ در حوضه او همدان- بهار. نسریه علوم و 
فنون کشاورزی و منابع طبیعی. علو مآب‌وحاک. سال چهاردهم. شماره ۵۳. صص ۹-۲۳. 

اعلمی. محمد تقی؛ عباسی. حبیبه و نیک سخن. محمد حسین؛ ۱۳۹۷. مقایسه دو روش متفاوت واسنجی و تحلیل عدم 
قطعیت مدل ٩۷۷۸‏ در برآورد میزان رواناب و بار مواد معلق حوضه صوفی چای. نشریه دان شآب‌وخحاک. جلد 
۸. شماره ۲ صص ۵۳-۱۶. 

اکبری مجدن حسین؛ بهره مند. عبدالرضا؛ نجفی نژاد. علی و واحد بردی. شبخ؛ ۱۳۹۲. شبیه‌سازی جریان روزانه 
رودخانه چهل چای استان گلستان با مدل 9۷۷۸ نشسریه پژوهش‌های حفاظت آب‌ونحاک. شماره ۳ صص 
۲۵۲-۹ 

انصاری. محمدرضا گرجی. منوچهر. صیاد. غلام عباس. شرفا؛ مهدی و حمادی کاظم؛ ۱۳۹۳. شبیه‌سازی رواناب 
حوضه آبخیز رود زرد با استفاده از مدل ۸ ۹۷۷. نشریه علوم و مهندس یآبیاری. جلد ۳۸. شماره ۶. زمستان 
صص ۹۷,۱۳۹۶-۱۰۷. 

بسالت پون علی اصغر؛ ایوبی. شمس ا..؛ حاج عباسی. محمدعلی و جلالیان. احمد؛ ۱۳۹۳. واسنجی و اعتبار سنجی 
مدل ۹۷۷۸ با استفاده از الگوریتم پرندگان برای شبیه‌سازی رواناب و رسوب در یک حوضه آبخیز کوهستانی؛ 
نشریه مدیریت حاک و تولید پایدار, جلد چهارم. شماره چهارم. ۱۳۹۳. صص ۲۹۵-۳۱۲. 

بصیری, مرضیه؛ اسفندیار پور بروجنی» عیبسی؛ اخوان» سمیرا؛ حمیدپون محسن و کمالی؛ اردوان؛ ۱۳۹۳. تأثیر مقیاس 
نقشه خاک بر مدل ٩۷۷۸‏ در شبیه‌سازی روند جریان رواناب ماهانه در حوضه آبخیز تویسرکان؛ نشریه مدیریت 
نحاک. تولید پایدار. جلد چهارم. شماره سوم. ۱۳۹۳. صص ۱۳۵-۱۵۳. 


سال هشتم برآورد میزان رواناب حوضه آبریز لنبران چای در استان ... قز 


بیات. مهراد؛ علیزاده حسین و مجردی, برات؛ ۱۳۹۷. جذب داده برای واسنجی - پیش‌بینی با اسستفاده از مدل 
شریه تحفیفات مناب مآب ایران. سال چهاردهم. شماره ۱. بهار ۱۳۹۷. صص ۰1-۱۲ 
حبیبی. علیرضا و گودرزی. مسعود؛ ۱۳۹۷. کاربرد مدل نیمه توزیعی ٩۷۷۸‏ در شبیه‌سازی رواناب حوضه حبله 
رود. نشریه علوم مهندسی وآبخحیز داری. سال دوازدهم. شماره 4۳. زمستان ۱۳۹۷. صص 1۰-4۹ 
رضازاده محمد سهیل بختیاری بهرام. عباسپو کریم. احمدی, مهدی و احمدی, محمد؛ ۰۱۳۹۲ شبیه‌سازی رواناب؛ 
رسوب و تبخیر- تعرق با استفاده از سناریوهای مدیریتی برای کاهش بار رسوب با استفاده از مدل .٩۲۷۸۲‏ 
نشریه علوم مهندس یآبخحیزداری ایران. سال دوازدهم. شماره 4۰. بهار ۱۳۹۷. صص 1۱-۵۰. 
زارع گاریزی آرش و طالبی. علی؛ ۱۳۹۵. شبیه‌سازی بیلان آب حوضه آبریز با استفاده از مدل 5۷۷۸ مطالعه 
موردی: حوضه قره سو استان گلستان. نشریه علوم مهندسی مناب عآب. سال نهم.۵۰ پاییز ۱۳۹۵. صص ۳۷-۵۰. 
شفیعی مطلق. خسرو؛ پر همت. جهانگیر؛ صدقی. حسین و حسینی. مجید؛ ۱۳۹۷. بررسی تغییر کاربری بر رواناب 
رودخانه مارون در ایستگاه ایدنک با استفاده از داده‌های سنجش ازدور و مدل ۸۵ 9۷۷؛ نشریه حفاظت منابع 
آب‌ونحاک. سال هفتم. شماره سوم. بهار ۰۱۳۹۷ صص ۷۱-۸۷ 
عارفی اصل, اکرم؛ نجفی نژاده علی و سلمان ماهینی» عبدالرسول؛ ۱۳۹۲. شبیه‌سازی رواناب و رسوب با استفاده از 
مدل ٩۷۷۸‏ در حوضه آبخیز چهل چای استان گلستان. نشریه مرتم وآبحیز داری. دوره ۰1۱ شماره ۳. پاییز 
۲ صص 1۳۲-۶۱ 
عینی. محمدرضا جوادی سامان و دلاور مجید؛ ۱۳۹۷. ارزیابی عملکرد داده‌های باز تحلیل شده پایگاه‌های اقلیمی 
جهانی [01) و 0۳58۴ ۱۳۳ در شبیه‌سازی هیدرولوژیکی مدل ٩۱۷۸۲‏ مطالعه موردی: حوضه آبریز 
مهارلو. نشریه تحفیقات مناب عمآب ایران. سال چهاردهم. شماره ۳ بهار ۱۳۹۷. صص ۳۲-4۶. 
غلامی, عباس؛ شاهدی, کاکا؛ حبیب نژاد روشن, محمود؛ وفاخواه مهدی و سلیمانی. کریم؛ ۱۳۹۲. ارزیابی کارایی مدل 
نیمه توزیعی ٩۷۷۸‏ در شبیه‌سازی جریان رودخانه‌ای (مطالعه موردی: حوضه آبخیز تالار استان مازندران)؛ 
نشریه تحفیقات آب‌وحاک ایران. دوره ۰4۸ شماره ۳ مهر و آبان ۱۳۹۲. صص 17۳-۶۷۱. 
کاویان. عطا الله. گلشن. محمد. روحانی. حامد و اسمعلی عوری. اباذر؛ ۱۳۹۲. شبیه‌سازی رواناب و بار رسوب حوضه 
آبخیز رودخانه هراز مازندران با بهره گیری از الگوی ٩۷۷۸۵۲‏ نشریه پژوهش‌های جغرافیای طبیعی. دوره 2۷. 
شماره ۲. تابستان ۱۳۹۶. صص ۰1۹۷-۲۱۱ 
کرمی. فریبا و بیاتی خطیبی مریم؛ ۱۳۹۸. مدل‌سازی فرس‌ایش خاک و اولویت بندی تولید رسوب در حوضه سد 
ستارخان اهر با استفاده از مدل‌های 91 ]1۷ و ۸ ۱۷۷. نشریه هیدروژئومورفولوژی. شماره ۱۸. سال ۵. بهار 
۸. صص ۰۱۱۵-۱۳۷ 
گودرزی. محمدرضا و فاتحی فر. آتیه؛ ۱۳۹۸. پهنه‌بندی خطر سیلاب در اثر تغییرات اقلیمی تحت سناریو 8.5 56۳ 
با استفاده از مدل هیدرولوژیکی ٩۷۷۸۵‏ در محیط 15 (حوضه آذرشهرچای). نشریه تحفیقات کاربردی علوم 


جغرافیایی. سال نوزدهم. شماره ۵۳. تابستان 4۸ صص ۹4-۱۱۷. 


۷ جغرافیا و مخاطرات محبطی شماره سی و یکم 


محمدی, مازیار؛ کاویان عطاا... و غلامی. لیلا؛ ۱۳۹۲. شبیه‌سازی دبی و نیترات در حوزه آبخیز تالار با استفاده از مدل 
نشریه پژوهشنامه مدیریت حوضهآبخیز. سال هشتم. شماره ۱۵. بهار و تابستان ۱۳۹۲. صص 1۵-1۰ 

مذن زاده روزبه؛ قهرمان, بیژن؛ ارشد. صالح و داوری قهرمان؛ ۱۳۹۵. بهبود مدل‌سازی حوضه آبریز از طریق تجمیع 
مولفه‌های اصلی هیدرولوژیک در مدل ۷۷۸ نشریه تحفیقات مناب عمآب ایران. سال دوازدهم. شماره ۳. پاییز 
۵ صص 1۵5-۷۹ 

ناصرآبادی فژاد؛ اسمعلی عوری, اباذر؛ اکبری» حسین و رستمیان رخساره؛ ۱۳۹. شبیه‌سازی جریان رودخانه با 
استفاده از مدل ٩۷۷۸‏ (مطالعه موردی: رودخانه قره سو اردبیل). نشریه پژوهشنامه مدیریت حوضه آبخحیز. سال 
هفتم. شماره ۱۳. بهار و تابستان ۱۳۹۵. صص ۵۰-۵۹. 


نظری پویاء هادی؛ کردوانی» پرویز و فرجی راد عبدالرضا؛ ۶ واسنجی و ارزیابی عملکرد مدل‌های هیدرولوژی 
1۳۱/۵۵۵ و 59۸ در شبیه‌سازی روانات. شریه تحلیل فضایی محاطرات محیطی. تحبال دوم. شماره ۲ 
تابستان ۱۳۹۶. صص ۹-۱۱۲. 
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0ص حمتاهطناه ۱۷۱۵01 :015مصظ 0مناماتتافیل ۵۶ متتعبظ م۲ ,1992 ,رعاصاظ به 220 بیک بقه6۷ظ 
279-۰ :6۰ ۴۳۵۵۰ 2۲۵ .ط0ا00عن همه 

۵۵ و ۵001و 9۷۷۸۵۲ فطا ۶ صمتاجناه۳۷ ,2004 و.تممصطصمطمصصتطاو به قصع ۱۷۷۰ .1 بتاطل) 
۸٩۸۵,‏ ظ) 0۶ فهمام1:۵282 .0صه وه۱۷ ۵۶ ممنعع مومع مردام ععمصطلعنم عطا مد اصعصمم‌مومع 
1057-0۰ :(47)4 

۵۴ 2۲66 2۵0 0106 ۱۷۲۵0۵111882 .2009 واط ۷202 یک منطو ناک عتامجرود‌ماطام ,۱۷ ۳۵۲۵۲۱۵۲21 
486-1۰ :23 ۲۲۵۵۵8569 ۲۱۲۲0۱021621 .موی ما اتازصاهانه 27 وععتتا0وع۲ 

٩۳/2] ۴0001 ۷716۶ 0‏ عطا عصتادع16 .2017 ,.[ ۲۵۲۵2 م2مطمصهو رو عمه* ناه ب۲ آتامصظ ,۷۲ ط0فونای) 
٩0۵0121 ۲۵5010610085. ۷۲۷۸۲۵۲, 901(4‏ امعم 0۴ 0262 تمطاهع ۷۷ 

۵۲ 01۲01021021 ۱ مط) عصتاع۱۷۵0 .2007 وی رعصمل من ره رل ,ان صمصصفط :2۵۵,۱۷۲ و۷ .> ولآ 
258-۰ :337 .۳۱۷۲0 0۴ .1 بقعظ ۷۷۵۵۲ ط1 ععصفطام0 ص12 ۵۶ 

,4132-0۰ ,(2-4) 333 ۱۱۷۲۵۱02۷ ۵۶ امطرا0ل 5۱۷۷۸۵۲۰ عملوه 60طوماه 20 رعما0ل بط عمتا[م۱۷]۱00 

61 .2007 ,با ۷۵100 ,بلط امصصته۲۱ بلط تمصعصاظ ,۱۷۲۷۷ 167 م2 ۷ ,10 10مصنه رل ۱۷۱۵:1281 

۰ ۱۷/۵۲۵۲۹۲۵0 1 پرمهتاهعع ۵ همهم ممتان ماود ت1۵ فممتامل‌تیام مماهتااه ام 

885-00۰ :(5063 ,تاظضففه مطا ۵۶ عومتاممعصم] ]1 

0 2 0۴ 211020107 ۴۷ .1999 ربق رعفطلم۴ 2240 ریک ر‌طه1ععظ۲ رما 00056 ,.۷۰0 1۷۲0۲۵۷1۱۵۷ 

96 ,۱۷۲۵۲۵۵۲0108 ۳۵۲۵5۲ 4صنه اوسانه نهد ممماججووهان اهلومرگ )عصتجعه ۱۷۲۵۵1 170102161 ۲۱ 

99, 257-7. 

۷۷۸۲ 001060 ۸ .2008 بیظ ماله0۷ع19نک مه ۳.۷۷۲۰ ,۱۷۲210 مب۷.ظ مهطاصهمل 

صهصارن .هم فممعوما :صتففظ ون تصهع‌مصهط ما مملام۱۷]۵0 ۷۱۵۱0 )ماع مه,ممتاهعزامج خر 

33-2 00۰ 2۰ ۱۱0۰ 2۰ ۷۵۰ .عصمعصنعصظ آماممجصصه110 ۱۳۷ 

0 01۱ .2004 .۲,۷۷ ,عصنک 20 بل مقفتصقهلال1 ۷۷ ویک روتنک ریت6 بلامصظ مساگه بطه‌کاله 
نماض .عمر16. .2005 مملوزع۷ 1۳0۵۵۵0 ۰۲۳060۵6621 -۰1001 )معصرووع2۱۵۲۰۸8۹ ۱۷۷ 
4 ۰ 6۲۷166 1۵562101 


سال هشتم برآورد میزان رواناب حوضه آبریز لنبران چای در استان ... ۷۹ 


1 .2010 ,,۷۰۵ ,قا2اص 1و1 20 یه ب6101)21) ریما .۲) مفهکللت) وت) ,۲۵2125 .۷ ,۲16101725 
1 26۲ ۷۷ 20 9011 معط عصلوها متعهط 912760 مضه ج ص1 همتام0مو «اتلمتان م۳2 2۳0 
٩0۰ 1۰ 00۰ 2741-6۰‏ .250 ۷۵۱۰ .12۵5۵111321108 ۲۰( ۸ 5۷۷) 1001 

ره ۷۷ ,1001۵ظاظ رصتفحظ ۷۵۲ ۱زا مباظ بعمهن] عط مد ممل)م۲۸1۵11012 ۷۷۵۵۵۲۹۴۵0 ۲۵۲ ط010)010عترظ 
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۰ 10601۲۵۱6۱6 1 ۵۵۵ ۲0 2 0 مم۵۵0۵ ۵۷ 0۲60102 ۶0۶ 00061 ۸ .1972 .1.1 160816[ 
.1204-1۰ :8 2700عوعک وعععتا0عمع ما2 ۱۷۷ 

صمصصنلعی ۵ عصتام۱۷۵0 .2010 ,. بخ رمقعمام]۷ 200 بیع رصقو۷2 صاتق ر.ظ ملطهع هن بت .9 مطعماه 
۶ 00/۵1 1۸۷۷۲۴۸ ,۱۵06۱ ۷۷۲۸۵۲ عملوتا ۱1001۵ ,لعصوع۲ه ۷۷ 0ععناهن)-لموزمض 1۳0 ۱۷1۵10 
514-۰ ,(46)3 ,۸9001200 عناوم ماع ۷۷ ج۲162همصخ مطع 

0 ۵01060 ۵۴ همتام۰۱۷۲۵0 .2010 وی ,۱۷۵1۵556 2040 یک رصهوه۷ 1101و ر.ظ بتطهع هن یوت مفلامصتطا 
۵۲ ۸۱۱۵۲۱۵۵0 ۰1۵ .۱۷۲۵06۱۰ ۶۷۷۸۲ م۱۲۹1 ۱1۳۱۱۵۳۲1۸ ,4عطوته2 ۷۷ 0عوناهت-نمعزدظ ۶۲۵ 
314-5۰ .00 .3 ۱۲0۰ .46 ,۷۵۱ .۸99001210108 قعععا0وم 

۷ صلوهها هعتهعهیان ما موم ارمو ۵۶ ممتهتاله۰۳۷ .2007 و.ظا به911۷7 20 .۵,۲ روماصهو ی ب۷2 1۱ 
۰ 2 ۱۱0۰ .1 ۷۵۱۰ .عم1زععمزممظ آه)ممجصمن عرص موهازرنا ۵061۰ 0عقهه همتقفمهو مامصصع 20 وآت) 
.44-0 

20 ممهتاه وصتاملم]۷ .1983 وال امعل رظ.6) تفای ونان ۱۷۲۵۲۲۵ ب06.ظ طوما 
معاتاتاف! موم ومممتا0فمع ماج ۷۷ .عمطوماه۳ ۵۲۵9060 همتمماف-تاومماد صم اتمه 
۰ 017 مک ۵ واتونه۷ ۲۳ .142 ۱۵۵0۲۲ 

هصلوه ۵061 ۱۷۷۸۵۸۲ طا ممتاهاتصلووه 0242 س«ماصصهمو .2015 وی تاملتمو 0ص با تمافالظ بنل نا 
1-4 ۲۱۷۲0۱08۵7 0۶ ما0[ ۲۱۱۵۵۲۰ صهصصلهک 060هع))ر۴ 

عصلوه تقو ۷1116*ععمصصعع فطع ۵۶ عملاع0مهه 0عحورع2 ۷۷ .2004 ,.عملمصصو م9 یه 6 0صع ۵ .ظ مصموان 1[ 
,110۷ 0 ما60 عطا تم ممتاهلتاه۱ مه ممتاه‌طانلهه رکصممون0۵۷610 ۱۷۵061 :5۱۷۷۸2000 
۵01 ۵۵0۲ 160۳01621 ,۲۵9۵۲۷۵1۲ ۷1۱16ععمصصهت) فط مه ع0مرعصهتا دبارمطامومطم قمع عصمحطللعو 
٩. ۰‏ ب262ط10 و۷ اامصمن رهمتهمممتممص مممصممز نم مه 011 ۵ 


